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Abstract 

Starting from the relativistic gauge-invariant quark-antiquark Green function we 
obtain the relevant interaction in the one-body limit, which can be interpreted as 
the kernel of a non-perturbative Dirac equation. We study this kernel in different 
kinematic regions, reproducing, in particular, for heavy quark the potential case 
and sum rules results. We discuss the relevance of the result for heavy-light mesons 
and the relation with the phenomenological Dirac equations used up to now in the 
literature. 
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1 Introduction 



In the last years a lot of efforts has gone into the study of systems involv- 
ing at least one heavy quark (Q). On the experimental side there are a lot 
of data on the heavy mesons states (QQ), new data already measured [1] 
and a great expectation for the ones to come on the heavy-light states (qQ). 
On the theoretical side the situation is the following. The dynamics of the 
systems composed by two heavy quarks is quite well understood in terms of 
potential interaction (static and relativistic corrections) [2-8] obtained from 
the semirelativistic reduction of the QCD dynamics (for lattice studies see 
[9]). In the heavy-light case it turns out useful to take advantage of the heavy 
quark symmetries [10]. Heavy-quark symmetry implies that, in the limit where 
rriQ ^> Aqcd, the long-distance physics of several observables is encoded in 
few hadronic parameters, which in general can be defined in terms of operator 
matrix elements in heavy quark effective theory (HQET). For some recent re- 
views we refer the reader to [11]. In this framework a systematic expansion can 
be done in the small parameters Aqcd/^q and a s (mg). In the limit mg — > oo 
all the physics can be expressed in terms of a small number of form factors 
depending on the light quark and gluon dynamics only. Then, in HQET the 
heavy-light meson mass is given by 

m M = ran + A + O \ ] corrections (1) 

\ m Qj 



The parameter A represents contributions coming from all terms independent 
of the heavy-quark mass itlq; it is one of the non-perturbative parameters of 
the HQET, which have a similar status as the vacuum condensates in sum 
rules and QCD phenomenology. A can be fixed on the data, its actual calcu- 
lation however needs a dynamical input. Of course the dynamics of the light 
quark is inerehently non-perturbative. Some approaches resort to dynamical 
calculation via phenomenological potential models [12,13], sum rules [14] or 
relativistic phenomenological equations [15]. To have a well founded calcula- 
tion of A is of great importance since its value affects the determination of 
many phenomenological quantities (cf e. g. [16]). 

In this letter we address the question of calculating the non-recoil corrections 
to the heavy-light mesons (A) via a Dirac equation justified by the QCD 
dynamics. Our starting point is the quark-antiquark gauge-invariant Green 
function taken in the infinite mass limit of one particle. The only dynamical 
assumption is on the behaviour of the Wilson loop. The gauge invariance of 
the formalism guarantees that the relevant physical information are preserved 
at any step of our derivation. In this way we obtain a QCD justified fully 
relativistic interaction kernel for the quark in the infinite mass limit of the 
antiquark. This kernel reduces in some region of the physical parameters to 
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the heavy quark mass potential, and leads in some other region to the heavy 
quark sum rules results, providing in this way an unified description. In the 
light of our result we scrutinize the phenomenological Dirac equations used in 
the literature and give an answer to the old-standing problem of the Lorentz 
structure of the Dirac kernel for a confining interaction [15,17-19]. 

There are many possible applications of the obtained result, like the study of 
relativistic properties of the spectrum (as much relevant as the quark is light) 
and the calculation of heavy-light meson matrix elements and form factors 
(e. g. the Isgur-Wise function). Finally, this work can also be intended as a 
step forward both in the direction of a theory derived two-body relativistic 
interaction, both in the direction of a generalization of the sum rules approach 
with the inclusion of a finite gluon correlation length. 



2 The one-body interaction 



The quark-antiquark Green function is given in quenched approximation by 

G inv (x, u, y, v) = (Tr i SW(x, y; A)U(y, v) i S^(v, u; A)U(u, x)) , (2) 



where the points x,y,u,v are defined as in Fig. 1, ( ) means the norma- 
lized average over the gauge field A^, is the fermion propagator in the 
external field A^ associated with the particle % and the strings U (y, x) = 

P exp |ig J ds(y — x) fl A^(x + s(y — x))\ are needed in order to have gauge 

invariant initial and final bound states. A very convenient way to represent it is 
the so-called Feynman-Schwinger representation (see [20,21] and refs. therein), 
where the fermion propagators are expressed in terms of quantomechanical 
path integrals over the quark trajectories (zi(ti) and ^2(^2)) 
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i J dh^F^( Zl ) i J dt 2 9 -a$F^(z 2 ) ^ (2) 
xe o e o (_j ^ +m 2 )\. (3) 

From Eq. (3) it emerges quite manifestly that the entire dynamics of the 
system depends on the Wilson loop: 

W(T; A) = TyP exp |ig j> dz^A^(z) j , (4) 



being T the closed curve defined by the quark trajectories and the endpoint 
strings U(y,v) and U(u,x). 

Let us assume that the antiquark moving on the second fermion line becomes 
infinitely heavy. The only trajectory surviving in the path integral of Eq. (3) 
associated with the second particle is the static straight line propagating from 
v to u. The corresponding Wilson loop of the system is represented in Fig. 
1. As already noted in [22] in this case it turns out to be very convenient to 
choose the following gauge condition (sometimes called modified coordinate 
gauge): 

A^(x , 0) = 0, x>Aj{x Q , x) = 0. 

Thanks to this gauge choice it is possible to express the gauge field in terms 
of the field strength tensor, 

i 

A (x)=J da x k F k0 (x ,ax.), 
o 

Moreover the only non-vanishing contribution to the Wilson loop is given by 
the quark paths connecting x with y, and we have 

W(T; A) = TrP exp l^ig J dsM^z) J . (5) 



We stress that the choice of the gauge is in this approach really arbitrary and 
motivated only by convenience. Being the formalism completely gauge invari- 
ant, by handling properly we would obtain exactly the same results within 
any gauge. 
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x = (— T/2, x) 



(T/2,y) 



u=(-T/2,6) 



« = (T/2,0) 



Fig. 1. The Wilson loop in the static limit of the heavy quark. 



As showed in [20,21] in order to evaluate Eq. (3) we need to know the Wilson 
loop average over the gauge fields. We evaluate it via the cumulant expansion 
described in [4]. Keeping only bilocal cumulants we obtain: 



y y 



(W(F, A)) = exp - 9 - J dx'» J dy'^D^x', y') , 
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D lxu (x,y) = x k y l J daa n ^ J d(3 (3 n ^ (F k ,(x° , ax)F lu (y°, (3y)), 



(6) 



where n(0) = and n(i) = 1. Assumption (6) corresponds to the so-called 
stochastic vacuum model and has been very successful in the last years either 
in applications to potential models as well as in the study of soft high energy 
scattering problems (for some recent reviews see [23]). Inserting expression (6) 
in Eq. (3) and expanding the exponential we obtain the following expression 
for the propagator Sp of the quark (So is G; nv "projected" on the first fermion 
line; the second quark is irrelevant in the infinite mass limit, playing the role 
of an external source): 



S D = S + S KS + S KS KSo + 



(7) 



So is the free fermion propagator. Taking into account only the first planar 
graph (since we are interested only in contributions proportional to the gluon 
condensate), we have K(y',x') = , y u So(y',x')'-f' J 'D f j, u (x',y'). A graphical rep- 
resentation of K is given in Fig. 2. Eq. (7) can be written in closed form as 
Sd = So + SoKSd (or in terms of wave-function, {j> — m — iK)ip — 0; m = mi). 
Therefore, K can be interpreted as the interaction kernel of the Dirac equation 
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associated with the motion of a quark in the field generated by an infinitely 
heavy antiquark. 



X 



X 



y 



y 



p \ k ■ q 



D 



Fig. 2. The interaction kernel K . 



Assuming that the correlator (F /iA (a;)F lyp (t/)) depends only on the difference 
between the coordinates, we define: 

(F k ,(x°, ax)F lu (y , /3y)> = - y°, ax - (3y). 

With this assumption K can be written in momentum space as (see Fig. 2 for 
the definition of the momenta): 



K( q ,p) = -g 2 (2 7 v)5(p°-q ) 

+00 i 1 

x J dr J daa nM J dp -— — J d 3 r e^-^ r 

-00 

x Y {#(-T)A + (t) 7 e- i(p0 - £t)r - #(T)A_(t) 7 e- l(p0+£t)T } Y 
x/ fc ^(r,(a-/3)r), (8) 

where t = ({3p-aq)/(f3-a), E t = Vt 2 + m 2 and A±(t) = — 

Equation (8) is our basic expression. It contains the perturbative interaction 
up to order g 2 and the non-perturbative one carried by a single insertion of 
a second order cumulant. From now on we want to focus our attention only 
on the purely non-perturbative interaction. The Lorentz structure of the non- 
perturbative relevant part of f p \ up is 

/JWs) = He { -Q^(g»»gx P - g, P gxu)D(x 2 ), (9) 



"c 
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where (F 2 (0)) is the gluon condensate, ll c the identity matrix of SU(3) and 
D is a non-perturbative form factor normalized to unit at the origin. Lattice 
simulations have showed that D falls off exponentially (in Euclidean space- 
time) at long distances with a correlation length a -1 ~ (1 GeV) -1 [24]. This 
behaviour of D is sufficient to give confinement, at least in some kinematic 
regions [4]. Moreover we notice that, if f^ vp ~ e~ MT , then we have 

K ~ W k ~k l / rf3r el(P " q) '^^bo + MhXt(0, (a - P)r). 

The main effect of the finite correlation length a -1 seems to consist, therefore, 
in a shifting of the pole in the inserted free fermion propagator. 

In the following we will study expression (8) for different choices of the pa- 
rameters which are the correlation length a, the mass m, the binding energy 
(po — m) and the momentum transfer (p — q). 

A. Heavy quark potential case (m > a > p — m) 

If we assume a to be bigger than the binding energy (p — tti) and smaller 
than the mass m of the quark, since a ~ 1 GeV, the quark turns out to be 
sufficiently heavy to be considered non-relativistic. In order to obtain the 1/m 2 
potential we can neglect the "negative energy states" contributions to (8) by 
writing 



K(q,p)~-g 2 (27r)5(p -q°) 

+ OC 1 1 



x J drj daa nM J d/J/^M — — J d 3 re i(p - q) ' r 



x/A + (t) 7 Va(r,(«^)r). (10) 

Now, inserting Eq. (9) and by means of usual reduction techniques [25], we 
obtain up to order 1/m 2 the static and spin dependent potential 



^)^ 2 ^ / drJdX(r-X)D(r 2 -X 2 ) 

C -oo 



where we have make use of the change of variable A = (a — (3)\r\. This result 
agrees with the one body limit of the potential given in [4,6]. In particular for 

(F 2 (0)) + f° °f 

r — > oo identifying the string tension a = g 2 — / dr dX D(t 2 — A 2 ) we 

— JL _i V s* J J 



— oo 
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obtain the well-known Eichten and Feinberg result [2], 



(12) 



(F 2 (0)) r °r 

where C is the constant term g 2 / dr dXX D{r 2 — A 2 ). We observe 

V c J J 

-co 

that the Lorentz structure which gives origin to the negative sign in front of 
the spin-orbit potential in (12) is in our case not simply a scalar (K ~ ar). 
We will discuss this point in more detail in the conclusions. 

B. Sum rules case (a < p — m, a < m) 

Let us consider now the case in which the binding energy of the quark is bigger 
than the correlation length, which can be considered zero respect to all the 
scales of the problem. In the literature is usually referred to this case as the 
non potential case [4]. Since /^(x)^/^(0) = U c (F 2 (0))/24iV c (^ Ap - 

Q^gxu), we have 



K(q,p) ~ -g 2 (2ir)5(p° - gO)^ <^M(^ H - g^ 9uk ) 



c 



x J daa n ^ J d/J/3»M__ ( 7 ^ (p) 7 M (27r) 3 5 3 (p - q)) . (13) 
o o 

In particular from Eq. (13) we obtain the well-known leading contribution to 
the heavy quark condensate [26]: 

d 4 p f d 4 q rc 1 («F 2 (0)) 



««)-/w/w a{SlW,t( " )S * )} 12 ™ 



C. Light quark case (a > m) 

Since we have reproduced the known results concerning heavy quarks, Eq. 
(8) should maintain some physical meaning also by considering heavy-light 
mesons with a strange quark (like D s and B s ). In this case the light quark 
mass is smaller than a: m s ~ 200 MeV < 1 GeV. Actually the case a > m 
has to be considered as the only realistic one concerning heavy-light mesons. 
Under this condition either the exponent (po — E t ) as well as (po + E t ) can be 
neglected with respect to a. Therefore we have: 



K(q, P )~-g 2 (2ir)5( P °-q ) 
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+ 00 1 1 

x J dr J daa nM J d(3 (3 n(u) -^-^ J d 3 re l(p - q) ' r 



x 7' 7^(r, (« - 0)r). (14) 

We observe, as an appealing feature of this expression, that in the zero mass 
limit it gives a chirally symmetric interaction (while a purely scalar interaction 
breaks chiral symmetry at any mass scale). On the other hand, by considering 
only its static contribution (i. e. neglecting all momentum dependent terms), 
we obtain (using the same definition of the string tension given previously): 

K(q, P )\ static = 9 2 (2n)8(p° - q°) J dhe^'^-ar. 

The factor 5/3 in front of a which arises naturally in this approach under the 
considered physical conditions, seems to supply an explanation for the fact, 
observed by many authors [15,17,18], that in order to reproduce the D s and B s 
spectra from a Dirac equation with scalar confinement a string tension almost 
twice respect to the usual value is needed. 



3 Conclusions 



In the literature, also recently, a Dirac equation with scalar confining kernel 
(i. e. K ~ err) has been used in order to calculate non-recoil contributions to 
the heavy-light meson spectrum [15,17,18]. The main argument in favor of this 
type of kernel is the nature of the spin-orbit potential for heavy quarks. This 
turns out to have a long-range vanishing magnetic contribution (according to 
the Buchmuller picture of confinement) and is completely described by the 
Thomas precession term. This situation is compatible with a scalar confining 
kernel. However, assuming more sophisticated confinement models with a big- 
ger sensitivity to the intermediate distance region, the spin-orbit interaction 
has no more such a simple behaviour. In particular there show up non zero 
corrections to the magnetic spin orbit potential. Moreover, the velocity depen- 
dent sector of the potential seems not to be compatible with a scalar kernel 
(we refer the reader to [6] for an exhaustive discussion). Therefore also from 
the point of view of the potential theory there are strong indications that the 
Lorentz structure of the confining kernel should be more complicate that a 
simple scalar one. This emerges also in our approach. The kernel (8) follows 
simply from the assumption on the gauge fields dynamics given by Eq. (6). 
In principle all the graphs constructed by inserting non-perturbative gluon 
propagators on the quark fermion line should be taken into account. Since we 
are interested only in terms proportional to o (or (F 2 (0))), we keep only the 
first one. When performing the potential reduction of this kernel in the heavy 
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quark case (A) we obtain exactly the expected static and spin-dependent po- 
tentials. Therefore our conclusion is that there exists at least one non scalar 
kernel which reproduces for heavy quark not only the Eichten and Feinberg 
potentials in the long distances limit, but also the entire stochastic vacuum 
model spin dependent potential. Moreover when considering a, the inverse of 
the correlation length, small with respect to all the energy scales (case B), the 
kernel (8) gives back the leading heavy quark sum rules results. It is possible 
to try to extend the range of applicability of Eq. (8) to more realistic cases, 
like D s and B s mesons where the light quark mass is smaller than the cha- 
racteristic correlation length of the two point cumulant (case C) . The relevant 
part of the kernel is also in this case not a simply scalar one. In the static 
approximation this kernel seems, indeed, to be compatible with the existing 
phenomenology. We notice, however, that the situation is quite different from 
QED where in the Coulomb gauge a static interaction emerges without any 
approximation since the transverse part of the propagator of the exchanged 
photon (the relevant contribution to the binding) vanishes when the second 
fermion line is taken infinitely heavy. In our approach a purely static contri- 
bution never emerges, being relevant to the binding not an exchange graph 
between the quarks, but the interaction with the background vacuum fields. 
Therefore any static approximation in a real heavy-light system, for which the 
light quark is expected to be far from a static one, appears doubtful. 
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